OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Proteomic Characterization of Plasmid pLAI for 
Biodegradation of PolycycMc Aromatic Hydrocarbons in 
the Marine Bacterium, Novosphingobium 
pentaromativorans US6-1 

Sung Ho Yun\ Chi-Won Choi\ Sang-Yeop Lee\ Yeol Gyun Lee\ Joseph Kwon\ Sun Hee Leem^, Young 
Ho Chung\ Hyung-Yeel Kahng^, Sang Jin Kim'*, Kae Kyoung Kwon'*, Seung I! Kim^'^* 

1 Division of Life Science, Korea Basic Science Institute, Daejeon, Republic of Korea, 2 Department of Biology, Dong-A University, Busan, Republic of Korea, 3 Department 
of Environmental Education, Sunchon National University, Sunchon, Republic of Korea, 4 Korea Institute of Ocean Science & Technology, Ansan, Republic of Korea, 
5 Department of Bio-Analytical Science, University of Science and Technology (UST), Daejeon, Republic of Korea 



Abstract 

Novosphingobium pentaromativorans US6-1 is a halophilic marine bacterium able to degrade polycyclic aromatic 
hydrocarbons (PAHs). Genome sequence analysis revealed that the large plasmid pLAl present in N. pentaromativorans 
US6-1 consists of 199 ORFs and possess putative biodegradation genes that may be involved in PAH degradation. 1-DE/LC- 
MS/MS analysis of N. pentaromativorans US6-1 cultured in the presence of different PAHs and monocyclic aromatic 
hydrocarbons (IVlAHs) identified approximately 1,000 and 1,400 proteins, respectively. Up-regulated biodegradation 
enzymes, including those belonging to pLAl, were quantitatively compared. Among the PAHs, phenanthrene induced the 
strongest up-regulation of extradiol cleavage pathway enzymes such as ring-hydroxylating dioxygenase, putative biphenyl- 
2,3-diol 1,2-dioxygenase, and catechol 2,3-dioxygenase in pLAl. These enzymes lead the initial step of the lower catabolic 
pathway of aromatic hydrocarbons through the extradiol cleavage pathway and participate in the attack of PAH ring 
cleavage, respectively. However, N. pentaromativorans US6-1 cultured with p-hydroxybenzoate induced activation of 
another extradiol cleavage pathway, the protocatechuate 4,5-dioxygenase pathway, that originated from chromosomal 
genes. These results suggest that N. pentaromativorans US6-1 utilizes two different extradiol pathways and plasmid pLAl 
might play a key role in the biodegradation of PAH in N. pentaromativorans US6-1. 
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introduction 

Environmental contamination by polycyclic aromatic hydro- 
carbons (PAHs) is a serious problem due to their toxic, 
carcinogenic and recalcitrant properties [1], and hence their 
biodegradation is an important process crucial for bioremediation, 
and understanding the breakdown pathways is an important part 
of developing clean-up technology. High-resolution analytical 
chemistry for metabolomics and high-throughput sequencing for 
genomics are essential for resolving PAH biodegradation path- 
ways. Recentiy, high-throughput proteomic approaches and 
integrated omics technologies have become important tools in 
the discovery of related proteins and enzymes [2,3]. For example, 
metabolite analysis and proteogenomic methods were recentiy 
used to fuUy understand the biodegradation of pyrene in 
Mycobacterium vanbaalenii PYR- 1 [4-7] . Proteomic studies have also 
been conducted on the PAH-degrading Sphingomonas sp. CHY-1 
and Mycobacterium sp. KMS [8,9] . 



Novosphingobium pentaromativorans US6-1 is a Gram negative 
halophilic marine bacterium able to utilize several PAHs, 
including phenanthrene, pyrene, and benzo[a]pyrene, as sole 
carbon and energy sources [10]. Genome sequencing of jV. 
pentaromativorans US6-1 has been recently completed and the 
genome database is accessible from the public NCBI database 
[11]. jV. pentaromativorans US6-1 contains two plasmids, pLAl and 
pLA2. Large plasmid pLAl possesses clustered putative aromatic 
compound degradation genes. The purpose of this study was to 
elucidate the PAH biodegradation pathways active in jV. 
pentaromativorans US6-1. Proteomic analysis of jV. pentaromativorans 
US6- 1 cultured in the presence of different PAHs was performed 
to identify biodegradation-related proteins and revealed the 
induction of low-molecular-weight (LMW) aromatic-hydrocarbon 
degrading genes located in pLAl. Importantiy, this strain uses a 
plasmid-born extradiol cleavage pathway (catechol-2,3 dioxygen- 
ase pathway) for the degradation of both PAHs and MAHs. In this 
study, we report on the role of plasmid pLAl in PAH 
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biodegradation and the physiological characterization of jV. 
pentaromativorans US6-1 using proteomic approach. Genomic 
studies on the biodegradation plasmids from several stains have 
been conducted previously [12]. However, proteomic character- 
ization of these extrachromosomal genetic elements has yet to be 
performed. This study reports the proteomic analysis of a 
biodegradation plasmid. 

Materials and Methods 

Bacteria cultivation and sample preparation 

jV. pentaromativorans US6-1 was cultured according to a method 
described previously [10]. A starter culture of bacteria was 
prepared by growing cells in marine broth (MB) at 30°Cto an 
optical density at 600 nm (ODgoo) of 0.8. Bacteria were harvested 
asepticaUy and equal amounts of bacteria (culture of 500 ml) were 
added to fresh BushneU-Hass broth (incorporating 30 g NaCl/L) 
(BD, USA) containing phenanthrene (50 ppm), pyrene (50 ppm), 
benzo(a)pyrene (50 ppm), benzoate (50 ppm), and /)-hydroxy- 
benzoate (50 ppm). ODgoo values of all cultures were checked at 
every six hours until 48 hours using spectrophotometer (Beckman 
coulter proteome Lab Du 800, USA). Bacteria were harvested 
after 12 or 36 h before suspending in 20 mM Tris-HCl buffer 
(pH 8.0) and then disrupting twice in a French pressure cell (SLM 
AMINCO, Urbana, IL, USA) at 20,000 Ib/inl Supernatants 
(crude cell extracts) were collected by centrifugation (1 5,000 x^, 
45 min) and subjected to oxygenase activity assay and proteomic 
analysis. Protein concentrations were determined by the Bradford 
method [13]. Enzyme activity assay and proteome analysis was 
conducted on the basis of same protem amount. 

Dioxygenase activity assay 

The activities of catechol 1,2-dioxygenase (GDI, 2), catechol 2,3- 
dioxygenase (CD2,3), protocatechuate 3,4-dioxygenase (PGD3,4), 
and protocatechuate 4, 5 -dioxygenase (PCD4,5) were determined 
using a UV spectrometer (Beckman Goulter Proteome Lab 
DU800, USA), as reported previously [14]. The activities of 



GD1,2 and GD2,3 were assayed by monitoring increases in 
concentration of cis, cu-muconate at A26() and 2-hydroxymuconic 
semialdehyde at A375, respectively. Activities of PCD3,4 was 
determined by monitoring the increase in concentration of P- 
carboxymuconate at A290 (absorbance decreased as |3-carboxy- 
muconate concentrations increased) and PGD4,5 activity was 
measured through the increase in 2-hydroxy-4-carboxy-muconic 
semialdehyde at A410, respectively. For each assay, one unit (U) of 
enzyme activity is defined as the amount of enzyme producing 
1 mmol of product per min. Activity assay of each sample was 
conducted at least three times for each sample. 

One-dimensional gel electrophoresis (1-DE) and in-gel 
digestion 

SDS-PAGE and in-gel digestion was performed as reported 
previously [15]. After electrophoresis and Goomassie blue staining, 
ID-gels were divided into 10 fractions according to molecular 
weight. Sliced gels containing fractionated protein samples were 
digested with trypsin (Promega, Madison, WI, USA) for 16 h at 
37°G after reduction with 10 mM dithiothreitol (DTT) and 
alkylation of cysteines with 55 mM iodoacetamide. The digested 
peptides were recovered with extraction solution (50 mmol/L 
ammonium bicarbonate, 50% acetonitrUe, and 5% trifluoroacetic 
acid). The extracted tryptic peptides were dissolved in 0.5% 
trifluoroacetic acid prior to further fractionation by LG-MS/MS 
analysis. 

LC MS/MS analysis 

LG MS /MS analysis was performed according to a modified 
version of a previously pubhshed method [16]. Tryptic peptide 
samples were loaded onto a 2G-V/V trap column (Waters, USA) 
for the enrichment of peptides and the removal of chemical 
contaminants. Goncentrated tryptic peptides were eluted from the 
column and directed onto a 10 cm x75 \xm i.d. C18 reverse phase 
column (PROXEON, Denmark) at a flow rate of 300 nl/min. 
Peptides were eluted using a gradient of 0 — 65% acetonitrUe for 
80 min. AU MS and MS/MS spectra were acquired in a data- 
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Figure 1 . Induction of two extradiol dioxygenases (catechol 2,3-dioxygenase and protocatechuate 4,5-dioxygenase) in polycyclic 
aromatic hydrocarbon or monocyclic aromatic hydrocarbon cultivation. MB; Marine broth, BaP; Benzo(a)pyren, Pyr; Pyren, Phen; 
Phenanthrene, Ben; Benzoate, PHB; p-Hydroxybenzoate. 
doi:1 0.1 371/journal.pone.009081 2.g001 
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Figure 2. Cluster analysis of the proteomes of N. pentaromativorans US6-1 induced in polycyclic aromatic hydrocarbons. Proteome 
datasets were obtained from stain US6-1 cultured from four conditions (MB; Marine broth, BaP; Benzo(a)pyren, Pyr; Pyren, Phen; Phenanthrene) and 
were hierarchically clustered. Up-regulated or down-regulated proteins were indicated as red or green, respectively. Identified proteins were 
categorized into six groups (CI ~C6) according to differential induction in PAHs. Each group was divided by COG function (Table S2). Subgroup C ~ V 
in COG function were defined as follows; [C] Energy production and conversion, [D] Cell cycle control, cell division, chromosome partitioning, [E] 
Amino acid transport and metabolism, [F] Nucleotide transport and metabolism, [G] Carbohydrate transport and metabolism, [H] Coenzyme transport 
and metabolism, [I] Lipid transport and metabolism, [J] Translation, ribosomal structure and biogenesis, [K] Transcription, [L] Replication, 
recombination and repair, [M] Cell wall/membrane/envelope biogenesis, [N] Cell motility, [O] Posttranslational modification, protein turnover, 
chaperones, [P] Inorganic ion transport and metabolism, [Q] Secondary metabolites biosynthesis, transport and catabolism, [R] General function 
prediction only, [S] Function unknown, |T| Signal transduction mechanisms, [U] Intracellular trafficking, secretion, and vesicular transport, [V] Defense 
mechanisms. 

doi:10.1371/journal.pone.0090812.g002 



dependent mode using a LTQ;Velos ESI Ion Trap mass 
spectrometer (Thermo Scientific, Germany). Each full MS (m/z 
range of 300 to 2,000) scan was foUowed by three MS/MS scans of 
the most abundant precursor ions in the MS spectra. For protein 
identification, MS/MS spectra were analyzed by MASCOT 
(Matrix Science, http://www.matrix science.com). The genome 



sequence database of jV. pentaromativorans US6-1 (GI: 359402640) 
was downloaded from NCBI and used for protein identification. 
Mass tolerance of the parent or fragment ion was 0.8 Da. 
Cabamidomethylation of cysteine and oxidation of methionine 
were considered to be variable modifications of tryptic peptides in 
the MS/MS analysis. 
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Nano-UPLC-MS^ Tandem Mass Spectrometry and 
database search 

An alternative MS/MS analysis was conducted using a nano- 
ACQUITY Ultra Performance LC Chromatography™ 
equipped Synapt^'^ HDMS System (Waters Corporation, MA, 
USA) as described previously [17]. The flow-through peptides 
were applied to a nano-ACQUITY UPTC BEH300 C18 RP 
column (180 |J,m x250 mm, particle size, 1.7 ^m). Trypsin- 
digested peptide mixtures were loaded onto the enrichment 
column (180 |a.m x20 mm, particle size, 5 |J,m), which was 
equilibrated with mobile phase A (3% acetonitrUe in water with 
0.1% formic acid) to remove salts and concentrate the peptides. 
Flow-through peptides were directly applied to a nano- 
ACQUITY UPLC BEH300 C18 RP column (180 ^m 
x250 mm, particle size, 1.7 |J.m) at a flow rate of 300 nl/min. 
The step gradient was as follows: 3 — 40% mobile phase B (97% 
acetonitrUe in water with 0.1% formic acid) for 95 min, 40 — 70% 
mobile phase B for 20 min, and then a sharp increase to 80% 
mobile phase B for 10 min. MS data analysis was carried out 
with the continuum LC-MS^ data using ProteinLynx Global- 
Server (PLCS) version 2.3.3 (Waters Corporation). The criteria 
for protein identification used in the PLGS search engine were 
applied with a peptide tolerance of 100 ppm, a fragment 
tolerance of 0.2 Da, and a missed cleavage allowance of 1. 
Analysis of quantitative changes in protein abundance (>95% 
confidence based on peptide ion peak intensities observed in low 
coUision energy mode (MS) in a triplicate set) was conducted 
using Expression Software version 2(Waters Corporation). 

Cluster analysis of proteomic data and prediction of 
protein properties 

The emPAI values were used in the cluster analysis of all 
analyzed proteins, and the proteome dataset was z-transformed 
and median-centered normalized. Analysis of the proteome 
dataset was performed by Pearson correlation and average linkage 
hierarchical clustering by Cluster and Tree View [18]. All proteins 
identified by proteomic analysis were classified according to the 
cluster of orthologous groups (COG) functions, and their 
subcellular localization of the identified proteins was predicted 
by Cello (v. 2.0; http;//cello.life.nctu.edu.tw/) [19]. Prediction of 
trans-membrane topology was performed using Phobius (http:// 
phobius.sbc.su.se/) [20]. 

Purification and characterization of plasmid pLAl 

jV. pentaromativorans US6-1 was grown in LB broth containing 
10 g/L NaCl. Plasmid DNA was isolated using a standard alkaline 
lysis procedure, and purified on NucleoBond columns (Nucleo- 
Bond Plasmid BAC Maxi kit, Clontech, USA). Genomic DNA was 
shared by careful pipetting of the DNA solution up and down 
several times with a 200 \i\ Pipetman tip. Purified plasmid DNA 
was analyzed in a 1 % (w/v) agarose gel by clamped homogeneous 
electrical field (CHEF) gel electrophoresis (Bio-Rad, USA). CHEF 
parameters were set according to the manufacturer's protocol, 
including a circulating temperature of 14°C, electric current of 
350 niA, and a pulse time of 35 s for 28 h. After electrophoresis, 
the gel was stained with ethidium bromide solution for 30 min, 
washed with TBE for 30 min, and DNA bands were detected 
using the Gel-Doc System (Bio-Rad, USA). 
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Figure 3. Two extradiol biodegradation pathways of N. pentaromativorans US6-1. Two biodegradation pathways were identified by 
genomic analysis using KEGG program and the induced enzymes were confirmed by proteomic analysis of N. pentaromativorans US6-1. PAHs are 
assumed to converge into 1,2-dihydroxynapthalene, whereas, ring-hydroxylating dioxygenase and dihydrodiol dehydrogenase are induced in our 
PAHs culture conditions. Up: Catechol 2,3-dioxygenase pathway (large plasmid pLAl). Down: Protocatechuate 4,5-dioxygenase pathway 
(chromosome). 

doi:10.1371/journal.pone.0090812.g003 



Results 

Comparative analysis of dioxygenase enzymes in 
Novosphingobium pentaromativorans US6-1 in response 
to polycyclic and monocyclic aromatic hydrocarbons 

jV. pentaromativorans US6-1 were pre-cultured in MB to obtain 
sufficient biomass to determine dioxygenase activity and proteo- 
mic analysis. Approximately equal quantities of cells were 
transferred into PAH and MAH culture media. Although all 
culture has same cell mass, each has different ODgoo values 
because of different absorbance of PAHs and MAHs. OD values 
are as foUow; benzo(a)pyren (0.5146 at ODgoo), pyren (0.5612 at 
ODfioo), phenanthren (0.604 at ODgoo), benzoate (0.4242 at 
ODfioo), and p- hydroxybenzoate (0.4366 at ODgoo). After 12 h 



incubation, delta OD values of benzo(a)pyren, pyren, phenanth- 
ren, benzoate, and /(-hydroxybenzoate were +0.0101, +0.1461, + 
0.0628, -0.0347, -0.0199, respectively. After 36 h incubation, 
delta OD values of benzo(a)pyren, pyren, phenanthren, benzoate, 
and /)-hydroxybenzoate were +0.0013, +0.0254, -0.008, - 
0.0492, +0.5695, respectively. The bacteria were harvested after 
1 2 h and 36 h incubation for enzyme activity assay and proteome 
analysis. Highest delta OD values of benzo(a)pyren and pyren 
were +0.0101 and +0.1461 at 12 hours incubation, respectively. 
Phenanthren was continually increased until 48 hours (+0.0628). 
Unexpectedly, cell mass under benzoate culture condition was not 
increased, whereas degrading enzyme activities were significantly 
increased. The activities of four major dioxygenases (CD 1,2, 
CD2,3, PCD3,4 and PCD4,5) were assayed using protein extracts 
from jV. pentaromativorans US6- 1 cultured in the presence of three 
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Figure 4. Purification of Plasmid pLA1 of N. pentaromativorans 
US6-1. Open circled plasmid was indicated by an arrow. Lane 1; 
ProMega-Markers Lambda Ladders, lane 2; DNA size standards 
8—48 kb, lane 3; plasmid from W. pentaromativorans US6-1 cultured 
in LB broth, lane 4; plasmid from N. pentaromativorans US6-1 cultured in 
p-hydroxybenzoate media. 
doi:10.1371/journal.pone.0090812.g004 

PAH,s and two MAHs to determine which biodegradation 
pathways were induced. Analysis of the jV. pentaromativorans US6- 
1 genome indicated the presence of only extradiol oxygenase 
genes; however, many unspecified oxygenase genes were identi- 
fied. In attempting to determine if jV. pentaromativorans US6- 1 also 
expressed intradiol oxygenase activities, we selected the four 
dioxygenase enzymes (CD 1,2, CD2,3, PCD3,4 and PCD4,5) that 
we considered would cover most degradation pathways in aerobic 
cultures. No activity of the intradiol dioxygenases CD 1,2 and 
PCD3,4 was detected. Activity of CD2,3 was high in cells cultured 
in media containing phenanthrene (1.59 U/mg at 12 h cultiva- 
tion), benzoate (0.53 U/mg at 36 h cultivation), and ^-hydro- 
xybenzoate (0.31 U/mg at 36 hr cultivation), whereas, CD2,3 
activity was detected to some degree in cells cultured under all 
conditions (0.02-0.04 U/mg), including MB media (Fig. 1). 
Activity of PCD4,5 was only detected in cells cultured in p- 
hydroxybenzoate (0.20 U/mg at 36 hr cultivation) media (Fig. 1). 
This suggests that the CD2,3 pathway could be a primary pathway 
for the degradation of PAHs and benzoate, while /(-hydroxy- 
benzoate is broken-down via the PCD4,5 or CD2,3 pathway in jV. 
pentaromativorans US6- 1 . 

Proteomic analysis of Novosphingobium 
pentaromativorans US6-1 cultured with polycyclic 
aromatic hydrocarbons 

Cytosolic proteins were prepared for shotgun proteomics using 
1-DE/LC MS/MS from cells cultured in different PAHs. 
Approximately 650—718 proteins were identified in cells grown 
in the presence of phenanthrene, pyrene, and benzo(a)pyrene. 
Comparative analysis between the three PAHs was made with cells 
grown in MB as a control (Table SI and Fig. SI). Analysis 



revealed that 494 proteins were commonly induced by all culture 
media, with 26 — 32 unique to the different aromatic hydrocarbons 
used as sole carbon sources. The identified proteins were divided 
into six groups (CI— C6) according to cluster analysis with each 
protein group arranged according to COG functions (Fig 2). 
Enzymes of the PAH and MAH cataboUc pathways were clustered 
in group C2 and group C6, and were expressed at higher levels in 
the presence of phenanthrene. This was particularly noticeable for 
the degradation enzymes associated with 'secondary metabolite 
biosynthesis, transport and catabohsm [Q]' (Table S2). Proteins in 
group C6 were strongly induced in those cells grown in MB and 
the primary up-regulated COG protein group was 'Translation, 
ribosomal structure and biogenesis Q]'. The ribosomal proteins 
induced during growth in MB increased by more than 1.69-fold 
compared to PAHs. 'Proteins uncharacterized or unknown [R or 
S]' were relatively higher in C3 group proteins, which were 
abundant under benzo [a] pyrene culture conditions. 

A notable outcome of the proteomic analysis was the strong 
induction of genes originating from pLAl. This large plasmid 
encodes 199 genes, and approximately 27 are clustered, consid- 
ered to be a putative coding region for LMW aromatic- 
hydrocarbon degradation. These genes could regirlate biodegra- 
dation of bicyclic aromatic ring compounds through to tricarbox- 
ylic acid cycle metabolites, such as acetaldehyde. Approximately 
27 — 36 proteins coded by genes located on pLAl were induced by 
three PAHs. Among the putative biodegradation genes, 24, 19, 
and 12 proteins were up-regulated in the presence of phenan- 
threne, pyrene, and benzo(a)pyrene, respectively, compared to MB 
(Table 1). Comparative results showed that the greatest amount of 
protein induction occurred in the presence of phenanthrene, 
which is consistent with the results of the dioxygenase activity assay 
(Fig. 1). Enzymes catalyzing the breakdown of bicyclic aromatic 
compounds to acetyl CoA were strongly induced. Ring-hydrox- 
ylating dioxygenase, dihydrodiol dehydrogenase, putative biphe- 
nyl-2,3-diol-l,2-dioxygenase, and catechol 2,3-dixoygenase were 
included. Although the activity of catechol 2,3-dixoygenase in 
those cultures containing pyrene was low, proteomic results 
showed a strong induction of both this enzyme and others 
involved in related biodegradation steps (Table 1). The signal 
intensity for those enzymes produced by plasmid pLAl in the 
presence of benzo(a)pyrene was lower than that for the other 
PAHs phenanthrene and pyrene. 

These results were verified by proteomic analysis using Nano- 
UPLC MS, confirming the up-regulation of JV. pentaromativorans 
US6-1 biodegradation genes by phenanthrene, pyrene, and 
benzo(a)pyrene (Table 2). These results showed that the 
biodegradation genes located on pLAl play a major role in 
the utilization of PAHs as sole carbon sources in jV. 
pentaromativorans US6- 1 . 

Proteomic analysis of Novosphingobium 
pentaromativorans US6-1 cultured with monocyclic 
aromatic hydrocarbons 

In response to the differential induction of extradiol 
dioxygenases of JV. pentaromativorans US6-1 in the presence of 
benzoate and /i-hydroxybenzoate, a comparative proteome 
analysis was conducted. Approximately 1,475 proteins were 
identified. 126 were induced by benzoate, and 85 by p- 
hydroxybenzoate (Table SI and Fig. SI.). Genomic analysis 
revealed the presence of a CD2,3 gene and a ring-hydroxylating 
gene on pLAl, which are responsible for the initial step of the 
lower catabolic pathway of aromatic hydrocarbons. Notably, 
four PCD4,5 genes (small and large subunits) were found to be 
localized on the chromosome. Proteome analysis showed that 
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when jV. pentaromativorans US6-1 was cultured in /i-hydroxy- 
benzoate, two PCD4,5 and /)-hydroxybenzoate degradation 
enzymes were induced (Table 3 and Fig. 3). These genes (gene 
number NSU_3623-NSU_381 1) were clustered on contig 58 of 
the chromosome (Table S3). However, none of the three PAHs 
or benzoate induced the expression of PCD4,5 and p- 
hydroxybenzoate degradation enzymes, suggesting that the 
chromosomal-born protocatechuate pathway plays a role in 
the breakdown of /)-hydroxybenzoate, but not in the degrada- 
tion of PAHs or benzoate. Unexpectedly, the biodegradation 
genes on pLAl were expressed in response to j!)-hydroxybenzo- 
ate, despite having no direct role in the breakdown of this 
compound. 

Discussion 

The genomic sequence of the bacterium A' pentaromativorans 
US6- 1 , which is able to utilize PAHs as its sole carbon source, was 
reported recently [10,1 1]. However, the genes coding proteins that 
are important for the biodegradation of PAH have not been 
completely annotated, and until now, their function has remained 
largely speculative. Two putative clusters containing the genes 
necessary for the breakdown of aromatic compounds were 
detected in the biodegradation gene region of plasmid pLAl and 
contig 58 of the chromosome (Table S3). In this study, a 
proteomic analysis oi pentaromativorans US6-1 cultured in the 
presence of three different PAHs was conducted to examine the 
differential expression of biodegradative genes. The results 
indicate that PAHs strongly induced the expression of biodegra- 
dation genes located on plasmid pLAl (Fig 2 and Table 1), but 
that other biodegradation gene clusters on contig 58 were not 
induced. These results suggest that the biodegradation genes on 
plasmid pLAl are essential for the biodegradation of PAHs. Semi- 
quantitative proteomic analysis using emPAI revealed that tri- 
(phenanthrene), tetra- (pyrene), and penta-aromatic compounds 
(benzo(a)pyrene) induce a differential biodegradation capacity 
(Table 1 and 2). The total amount of biodegrading enzymes 
induced in the presence of phenanthrene was estimated to be nine 
times greater than that of benzo(a)pyrene (Table 1). However, we 
found a discrepancy between enzyme activity and the amount of 
CD2,3 expression induced in cells cultured in the presence of 
pyrene, and the reason for the instability of CD2,3 induced by 
pyrene-containing culture media remains unclear. Because genes 
that regulate the biodegradation of bicyclic aromatic compounds 
were identified on pLAl, genes for high-molecular-weight PAHs 
should be found on the chromosome. An analysis using the 
program pFam revealed that several putative aromatic compound- 
degrading genes are scattered throughout the chromosome 
(contigs 58, 54, and 55), although they were not induced 
significantly under our culture conditions. We also considered 
the possibility that the cytochrome P450 monooxygenases (CYPs) 
may provide alternative biodegradation pathways. Proteomic 
analysis showed that four cytochrome P450 proteins 
(NSU_2269, 2261, 2259, and 2257) were present, with the 
hypothetical protein NSU_2261 assigned as a cytochrome P450 
monooxygenases in uniprot Blast analysis. Induction level was very 
low, suggesting that CYP has no direct involvement with the 
biodegradation of PAHs (Table SI), and this was further 
supported by the genomic analysis of jV. pentaromativorans US6-1. 
The degradation of PAHs was catalyzed by enzymes with broad 
substrate specificities. Consequently, further biochemical assay 
and overexpression are required to determine related biodegrad- 
ing genes [21-23], coupled with a more accurate proteomic study 
and optimization of culture conditions to understand PAH 



metabolism in M. pentaromativorans US6-1. Plasmids similar to 
pLAl have been identified in Sphingomonas aromaticivorans F199 and 
Sphingomonas sp. strain KAl [24,25]. S. aromaticivorans F199 contains 
a large plasmid, pNLl, which possesses 186 ORFs and 79 genes 
that mediate catabolism or transport of aromatic compounds, such 
as mono-aromatic compounds (m-xylene and p-cresol) and bicyclic 
aromatic compounds (biphenyl and naphthalene). Plasmid 
pCAR3 from Sphingomonas sp. strain KAl contains a number of 
carbazole degradation genes. A comparative analysis of DNA 
sequence of pLAl and pNLl reveals that they are significandy 
homologous (more than 36%) and each of the ORFs in pLAl has 
more than 61% amino acid sequence homology with its 
corresponding ORFs (Fig S2). These plasmids are categorized 
into four groups according to the repA proteins, with pLAl from 
TV. pentaromativorans US6- 1 belonging to the Rep_3 superfamily. 

The proteomic characterization in this study revealed that jV. 
pentaromativorans US6-1 utilizes different pathways for the break- 
down of the two MAHs. Benzoate was degraded via the CD2,3 
pathway encoded by genes located on plasmid pLAl, whereas p- 
hydroxybenzoate was broken down by the PCD4,5 route. The 
lower induction of biodegradation enzymes coded by pLAl in the 
presence of MAHs such as benzoate and /i-hydroxybenzote 
suggests that their regulation was optimized by PAHs or their 
metabolites. This assumption should be confirmed by further 
studies. The genes on plasmid pLAl have been confirmed only by 
genomic sequencing and gene assembly. In a previous study, the 
presence of plasmid pLAl in jV. pentaromativorans US6-1 was not 
confirmed due to difiiculties in the purification of large plasmid (> 
200 kb). Here, we confirmed the presence of pLAl in JV. 
pentaromativorans US6-1 by CHEF gel electrophoresis, and identi- 
fied an open-circled plasmid of approximately >25() kb (Fig. 4). 
The presence of biodegradation genes on the plasmid was 
confirmed by PCR using the purified DNA as the template (data 
not shown). 

In conclusion, two extradiol degradation pathways mediated by 
genes located in the plasmid and chromosome were separately 
induced by PAHs and MAHs. The large plasmid pLAl plays a 
pivotal role in the degradation of bicylic aromatic compounds to 
acetyl CoA in jV. pentaromativorans US6-1. 

Supporting information 

Figure SI Total number of identified proteins of N. 
pentaromativorans US6-1 by LC-MS/MS analysis. A & B; 

The proteomes induced in PAHs and MAHs were identified using 
jV. pentaromativorans US6-1 genome database. Benzoate (Ben), p- 
hydroxybenzoate (PHB), phenanthrene (Phen), pyrene (Pyr), 
banzo(a)pyrene (BaP). C & D; the proteomes induced in PAHs 
and MAHs were identified using plasmid (pLAl) database of jV. 
pentaromativorans US6-1. 
(TIF) 

Figure S2 Comparative analysis of pLAl of N. pentar- 
omativorans US6-1 and pNLl of N. aromaticivorans 
pentaromativorans F199. Biodegradation genes were indicat- 
ed with red boxes. 
(TIF) 

Table SI Total identified proteins by LG-MS/MS 
analysis. 

PCLSX) 

Table S2 Clustering of total identified proteins of N. 
pentaromativorans US6-1. 

PCLSX) 
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Table S3 Putative biodegradation gene clusters of N. 
pentaromativorans US6-1 and their proteomic result 
according to PAHs and MAHs. 

PCLSX) 
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